Experimental study of the decay K\ — > n°vv 
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The first dedicated search for the rare neutral-kaon decay Kl — >• iy°vv has been carried out in the 
E391a experiment at the KEK 12-GeV proton synchrotron. The final upper limit of 2.6 x 10~ 8 at 
the 90% confidence level was set on the branching ratio for the decay. 

PACS numbers: 13.20.Eb, 11.30.Er, 12.15.Hh 



I. INTRODUCTION 

The study of rare kaon decays [l[ has played an im- 
portant role in the establishment of the standard model 
(SM) in particle physics. It has also been crucial for un- 
derstanding the phenomenon of CP violation Q. The 
rare decay K\ — > tt°vv [^-Q is a direct CP violation pro- 
cess caused by a flavor-changing neutral current (FCNC) 
with transition from a strange to down quark. 

The unique characteristic of the K\ — > tt q i/D decay 
is that the branching ratio can be calculated with very 
small theoretical uncertainties. In the SM based on 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix Q for 
quark flavor mixing, the K\\ — > ifivv branching ratio is 
predicted to be (2.49 ±0.39) x 10 -11 Q. The uncertainty 
of the prediction is dominated by the allowed range of 
the imaginary part of a CKM matrix element Vtd, which 
is determined by other measurements, and the intrinsic 
theoretical uncertainty is only 1-2% 

The determination of CKM matrix elements has been 
greatly improved in the past ten years by measuring var- 
ious B-decay properties Q; all measurements were con- 



sistent with each other within the standard model pa- 
rameterization. By measuring the K\\ — > iftvv decay 
precisely, we can check the consistency against the cur- 
rently predicted value in the SM. A deviation indicates 
new physics beyond the SM because of the very small 
uncertainty in deriving the imaginary part of Vtd from 
the branching ratio of the K\\ — > tt°i/D decay. 

The decay A'£ — > ifivv is one of the processes expected 
to have a significant impact on new physics searches be- 
cause it is an FCNC process that can proceed through 
loop diagrams, including the interactions at short dis- 
tance and large mass scales. Grossman and Nir Q 
pointed out the importance of studying new physics by 
using both the K\ — > ifivv and K + — > -k + vv decays, 
where the decay K + — > ti + vD is another rare kaon de- 
cay with small theoretical uncertainties [1,0, El- Various 
new physics models have been developed and used to pre- 
dict the AT? — > ii a vv and K + — > tt + vv branching ratios 

ana 

The best upper limit obtained in previous experiments 
was BR(K° L -> tt°vv) < 5.9 x 10~ 7 (90% CL) 0. The 
ultimate goal of our experimental study is to determine 
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the Kl — > ir°vD branching ratio with an accuracy less 
than 10% of the value predicted in the SM. The goal can 
be achieved by performing a series of experiments with 
improved and refined detection methods. 

The E391a experiment, which was carried out at the 
KEK 12-GeV proton synchrotron (KEK-PS), is the first 
step of this approach. The main objectives of E391a were 
not only to investigate the decay with a dedicated detec- 
tor, but also to test and confirm our basic experimental 
methods for K® L — > tt^vv at the highest possible sensitiv- 
ity. Data collection of the E391a experiment started in 
February 2004 and continued until November 2005, which 
was one month before the shutdown of the KEK-PS. The 
total running time was about twelve months and it was 
divided into three periods (Run-1, Run-2, and Run-3). 
Early results from the first and second periods were re- 
ported in Ref. [l3| and Ref. [HI, respectively. In this ar- 
ticle, we report the final results on the Kl — » ir°vv decay 
as obtained from E391a. 



II. EXPERIMENTAL METHOD 
A. Basic method of K { [ —¥ ir a vh> detection 

The key signature of the K\ — > ir°vv decay is detection 
of exactly two photons and nothing else. The tt° decays 
dominantly into two photons (77), and the neutrinos are 
undetectable. Because detection of the incident K\ is 
difficult, measurements can only be made of the energy 
and position of the two photons in a calorimeter located 
downstream of the decay region, without having di- 
rect information of the incident particle. Kinematic con- 
straints are weak for definitive identification of the de- 
cay. Instead, the decay must be isolated by eliminating 
all possible backgrounds. 

Because the signal mode is identified as the final 
state of two photons and nothing else, processes that 
make two or more photons can cause background events. 
The K\ decays such as K° L — > 7r°7r° and K° L 
7T + 7r~7r° become backgrounds when extra photons or 
charged particles escape detection. The K\ — > 77 decay 
can be a background source because it has only two pho- 
tons in the final state, although it is well suppressed by 
kincmatical constraints. Another background process is 
hadronic interactions of beam neutrons with the residual 
gas in the beamline or in detectors near the beam. Any 
7T° or rj produced in these interactions decay into two 
photons and produce backgrounds. Hyperons produced 
at the target can cause backgrounds through processes 
such as A — » 7r°n; these hyperon events are strongly sup- 
pressed because most of them decay in a 10 m long neu- 
tral beamline. 

The most important tool for reducing the background 
is a hermetic detector system to detect and veto extra 
particles. Because all of the other K\ decays, except 
for K*l — > 77 , are accompanied with at least two addi- 
tional photons or charged particles, the detector should 



be highly sensitive to photons and charged particles. 

The signature of K\ — > tt°vv can be provided by us- 
ing a small-diameter K® L beam (called a "pencil beam") 
and measuring precisely the energy and position of the 
two photons. Although flux is reduced with a pencil 
beam, it has several advantages. First, the beam hole 
at the center of the calorimeter, which compromises her- 
meticity, can be minimized. Secondly, the if£ decay ver- 
tex position (Zyrx) can be assumed to be on the beam 
axis. The vertex position, which is same as the ir° decay 
position due to the short lifetime of 7r°, is obtained from 
the kinematics of the two photons from its decay (See 
Sec. III.B). The transverse momentum of the ir° (Pt) 
with respect to the Kl beam axis is also obtained. 

The signal region for Kl — > ifivv decay can be de- 
fined with Zvtx and Pt- Requiring a sufficiently large 
missing transverse momentum eliminates contamination 
from the K\ — > 77 decay and also reduces the contami- 
nation of other K\ decays. Most of the K®'s decay into 
multiple particles that have low momenta in the Kl rest 
frame, and hence low Pt in the laboratory frame. The 
Zvtx should be in the region away from beam counters 
and P T should be in the range from 120 to 240 MeV/c, 
where the maximum momentum of 7r°'s in the A"£-rest 
frame is 231 MeV/c for K° L -)• ir^vv. 

Production of ir°'s through beam-gas interactions in 
the decay fiducial region are reduced by having a high 
vacuum. Hit-rates of the beam halo (mostly neutrons) 
with the surrounding detectors are minimized by a sharp 
collimation of the beam. In addition, as discussed in Sec. 
III. A and in the subsequent analysis of the data, mis- 
combinations of two photons from background processes 
and/or mis- measurements of energies and positions were 
found to be major sources of backgrounds in the E391a 
experiment. The photons were often produced from in- 
teractions of the beam halo and can lead to incorrect 
determinations of Zvtx- Some backgrounds from beam 
interactions can be reduced by detecting low-energy de- 
posits from recoil particles emitted after the interactions. 
More detailed descriptions of the detection method have 
been reported in the original proposal of E391a [lj|. 

Inefficiencies of photon and charged particle detection 
provide backgrounds. The decay K a L — > 7r°7r° is the most 
serious background source arising from an inefficiency be- 
cause it has only two extra photons in the final state. 
The inefficiency of photon detection was investigated in 
a series of experiments fl5| , which showed that the ineffi- 
ciency monotonically decreased with the incident photon 
energy. It was also found that photon detection with a 
very low energy threshold is necessary to achieve a small 
overall inefficiency, even for high-energy photons. If an 
extremely low energy threshold around 1 MeV was set, 
the backgrounds from other K^ decays were reduced to a 
negligible level within the sensitivity of the experiment. 

The following sections describe the application of the 
basic methods in the E391a experiment. 
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B. Beamline 

The production target of the K\ beam was a platinum 
rod with a length of 60 mm and a diameter of 8 mm. The 
profile of the primary beam was a = 3.3 mm and a 
i .1 mm along the X (horizontal) and Y (vertical) axes. 
The neutral beam was extracted at an angle of 4° with 
respect to the primary proton beam. The target rod and 
beamline elements were aligned along a straight line in 
the direction of the neutral beam. The total length of the 
beamline was 10 m. The neutral beamline consisted of 
two dipole magnets (Dl, D2) to sweep charged particles 
out of the beam, with field strengths of 4 T-m and 3 T-m, 
respectively, and six collimators (CI - C6) to collimatc 
the beam, as shown in Fig. [1] 
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FIG. 1. Schematic view of the neutral beamline: (top) the ar- 
rangements of the components and (bottom) the collimation 
scheme. C1-C3, C5, and C6 were tungsten collimators, assem- 
bled as a stack of cylindrical blocks 5-cm thick with circular 
holes of different diameters. Each collimator approximated 
the cone-shaped aperture as indicated by the A-E lines in the 
bottom figure. The Z = coordinate is at the center of the 
production target. 

Five collimators (C1-C3, C5, and C6) were made of 
tungsten. The C2 and C3 collimators in the upstream 
end were used to define the beam with the designed aper- 
tures, which were arranged to form a half cone angle 
of 2 mrad from the target center. The last two colli- 
mators, C5 and C6, were used to trim the beam halo. 
The most upstream collimator CI reduced the size of 
the beam immediately after the target without produc- 
ing a large penumbra. The total thickness of these five 



collimators was approximately 6 m. A thermal-neutron 
absorber, which was made of polyethylene terephthalate 
(PET) sheets containing Gd203 40% in weight, was used 
for C4. The aperture of C4 was set to be larger than that 
of the other collimators. 

Movable absorbers, made of lead (Pb) and beryllium 
(Be), were placed between CI and C2 to reduce the num- 
ber of photons and neutrons relative to the A"°'s. The 
absorbers were 10-mm-diameter rods with the lengths of 
5 cm and 30 cm for Pb and Be, respectively. The down- 
stream region, starting with a stainless steel window 100 
fim thick at the upstream end of C4, was evacuated to 
approximately 1 Pa. 

The primary beam on the target was monitored by 
a secondary-emission chamber (SEC) placed upstream 
of the target, and a target monitor (TM) that was 
a counter telescope which viewed the target center at 
90 degrees. The primary beam position at the target 
was adjusted with steering magnets by monitoring the 
K\ — > 7r 7r 7r° decay events. In the 37r° events, six 
photons were detected by a Csl calorimeter which will 
be described in Sec. II. C. The center of energy was de- 
fined to be r c = J^EiVi, where Ei are the photon en- 
ergies and Ti are the photon hit positions at the Csl 
calorimeter. Because the six-photon events were mostly 
K\ — > 7r°7r°7r decays, the peak position of the center of 
energy should be on the axis of the beamline. The po- 
sition of the neutral beam was maintained to be within 
0.2 mm from the center throughout the entire running 
period. The beam size at the Csl calorimeter, was also 
monitored by the distribution of the center of energy. Its 
diameter was a = 40 mm, which was consistent with the 
beam divergence. 

The peak momentum of the as determined from 
beamline simulations was 2 GeV/c at the exit of C6, as 
shown in Fig. [2] 

The initial Pt spread due to the beam divergence of 
2 mrad was approximately 4 MeV/c. The neutron-to- 
Kl ratio was 60 and the halo-to-core ratio was approxi- 
mately 10 -5 for both neutrons and photons with energies 
above 1 MeV, as shown in Fig. [3] By inserting a Be ab- 
sorber, the neutron-to-AT° ratio was reduced to 40, with 
roughly a 45% loss in K\ flux. Profiles of the K\ beam 
are shown in Fig. 

Punch-through muons were emitted in the direction 
parallel to the beam axis. Their position distribution 
was almost flat and the flux density was larger than the 
cosmic ray flux by roughly one order of magnitude. De- 
tails of the beamline has been reported elsewhere [13] • 



C. Detectors 

The E391a detection system was located at the end 
of the beamline. The detector subsystems were cylindri- 
cally arranged around the beam axis, and most of them 
were placed inside a large vacuum vessel, as shown in 
Fig. [5] From here on, the origin of the coordinate system 
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FIG. 2. Momentum spectrum of K\ at the exit of C6, ob- 
tained from beamline simulations. 



FIG. 4. Beam profile of K\ at the exit of C6 collimator, 
obtained from beamline simulations. Filled and open circles 
show the case with and without the Be absorber, respectively. 
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FIG. 3. Beam profiles of neutrons and photons above 1 MeV 
at the exit of C6 collimator, obtained from beamline simula- 
tions. 

is defined to be at the upstream end of the E391a detec- 
tor, as shown in Fig. [5j This position was approximately 
12 m from the production target. 

1. Csl calorimeter 

The energy and hit position of photons, such as the 
two photons that would come from the ir° decay in the 
Kl —> n°vv signal, were measured by using a calorime- 
ter placed at the downstream end of the decay region. 



As shown in Fig. [6l the calorimeter was made of 496 un- 
doped Csl crystals each with a dimension of 7 x 7 x 
30 cm 3 (Main Csl) [lH and assembled in a cylindrical 
shape with an outer diameter of 1.9 m. There was a 12 
cm x 12 cm beam hole at the center of the calorimeter. 
The beam hole was surrounded by a tungsten-scintillator 
collar counter (CC03) and 24 Csl modules with dimen- 
sions of 5 x 5 x 50 cm 3 (KTeV Csl) [Hj]. The gaps 
at the periphery of the cylinder were filled with 56 Csl 
modules with trapezoidal shape and 24 modules of lead- 
scintillator sandwich counters (Lcad/Scintillator Sand- 
wich). These modules were tightly stacked such that 
the gaps between them were less than 0.1 mm. Al- 
though the CC03 and the sandwich counters were used 
only in the veto system, the other parts, which con- 
sisted of Csl crystals, were used for photon measure- 
ments. The average light yield of the main Csl modules 
was 16 pe/MeV, where pe/MeV is the unit of the num- 
ber of photo-electrons emitted from the photomultiplicr 
tube (PMT) cathode for an energy deposit of 1 MeV in 
the detector. Details of the Csl calorimeter have been 
reported in Ref. [2(j . 



2. Charged veto counter 

A set of plastic scintillation counters (named CV) were 
placed in front of the Csl calorimeter to identify events 
that included the emission of charged particles, such as 
K\ — > -K~e + v decay [2l|. A total of 32 sector-shaped 
modules (outer CV) were placed at a distance of 50 cm 
from the front face of the Csl calorimeter. The outer CV 
modules were arranged to have overlaps between adjacent 
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FIG. 5. Detection system. 
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FIG. 6. Assembly of detector subsystems in the downstream 
end-cap. 



modules. The beam region from the outer CV to the Csl 
was covered by the inner CV, which was a square pipe 
formed with four plastic scintillator plates. The inner 
and outer CVs were closely connected with aluminum 
fixtures. In order to eliminate gaps between the outer and 
inner modules, the inner modules were extended to cover 
the edge of the outer modules. The edges of the outer 
CV and the inner CV were located close to the beam, 
and neutrons in the beam halo frequently interacted with 
them. 



3. Barrel counters 

The decay region was surrounded with two large lead- 
scintillator sandwich counters: the main barrel (MB) and 
the front barrel (FB), composed of 32 and 16 modules, 
respectively. The modules were tightly assembled with 
small gaps [22[ , as shown in Fig. [7J 

The lamination was parallel to the beam axis, and the 




FIG. 7. Cross-sectional end view of the MB (outer ring) and 
FB (inner ring). 



lengths of the MB and FB modules were 5.5 m and 2.75 
m, respectively. We developed a new type of plastic scin- 
tillator made of a resin mechanically strengthened with a 
mixture of styrene and methacrylate [23J. The scintilla- 
tor counters were fabricated by extrusion. The thickness 
of the MB and FB modules were 13.5Xo and 17.2Xo, re- 
spectively. The emitted light was transmitted through 
wavelength-shifting (WLS) fibers, which were glued to 
each scintillator plate with a pitch of 10 mm. The fibers 
had a doublc-cladded structure with Yll as a dopant. 
The green light emitted from Yll was read by a newly 
developed photomultiplier tube with a high quantum ef- 
ficiency for green light (EGP-PMT) 0]. For readout, 
each module was divided into an outer and an inner parts. 
The MB module was viewed from both ends, and the FB 
module was viewed from the upstream end. Figure [5] 
shows the light yields for 4 readout channels of one MB 
module, obtained from a cosmic-ray test. The light yield 
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and attenuation of the FB module were similar to those 
shown in Fig. [5] for the MB module. Details of the barrel 
counters has been reported in Ref. [25|. 
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FIG. 8. Light yield and attenuation for a MB module. 

A layer of plastic scintillator was placed inside of the 
MB modules to identify charged particles. It consisted 
of 32 modules and was called the barrel charged veto 
(BCV). A module was made of a 10-mm-thick plastic 
scintillator, and the light signal was read from both ends 
by an EGP-PMT through WLS fibers glued to the scin- 
tillator with a 5-mm pitch. 

4- Counters close to the beam 

Multiple collar-shaped counters, CC02 - CC07, were 
placed along the beam axis. Counter CC02 at the en- 
trance of the decay region was a lead-scintillator sand- 
wich of the Shashlik type, with WLS fibers piercing the 
lead and scintillator layers. Counter CC03 at the end 
of the decay region was a tungsten-scintillator sandwich. 
Counters CC04 - CC07 covered the solid-angle in the 
downstream direction. Counters CC04 and CC05 were 
lcad-scintillator sandwiches with WLS fibers glued on 
each scintillator plate at a pitch of 10 mm, which was 
the same as for the MB and FB. Counters CC06 and 
CC07 were made of SF5 lead-glass. The direction of lam- 
ination was perpendicular to the beam for CC02, CC04, 
and CC05, and was parallel to the beam for CC03. The 
signals from CC02 - CC05 were read by EGP-PMTs. 

Beginning with Run- 2, another collar counter called 
CC00 was installed in front of the FB and outside the 
vacuum vessel to reduce the effects of halo neutrons. It 
consisted of 11 layers of a 5-mm-thick plastic scintilla- 
tor interleaved with 10 layers of 20-mm-thick tungsten. 



However, CC00 did not significantly reduce the neutron 
halo because it was not placed as close to the beam as 
the other collar counters. 

A beam- plug counter, back-anti (BA), was placed at 
the end of detector system along the beam axis. In the 
first two data-taking periods, Run-1 and Run-2, the BA 
consisted of six superlayers, where each superlayer had 
six plastic-scintillator layers interleaved with lead sheets, 
and a single layer of quartz. In Run-3, the lead and 
plastic layers were replaced with PWO crystals as shown 
in Fig IH1 with the intention to separate electromagnetic 
showers from the neutron hits. 

A thin layer of plastic scintillator, the beam hole 
charged veto (BHCV), was placed in front of the BA. 
It was effective in removing the lf£ — > 7r°7r + 7r~ decay. 
The thickness of the BHCV was 1 mm for Run-1, and 
3 mm for Run-2 and Run-3. 



lead/scinti. module 
Quartz 



beam 




6x lead (2mm) + 7x scint. (5mm) 



PWO Quartz 




FIG. 9. BA used in Run-2 (top) and Run-3 (bottom). 

The upstream counters, CC02 - CC04, were placed 
inside the vacuum vessel, and the downstream CC05 - 
BA were placed outside the vessel. A vacuum duct, which 
was directly connected to the vacuum vessel, penetrated 
CC05 - CC07 and the vacuum region was extended up 
to the front face of BHCV. 

The size and basic parameters of CC02 - CC07, 
BHCV, and BA are summarized in Tabic |TJ 



5. Vacuum 

The vacuum pressure in the decay region had to 
be maintained below 10 -5 Pa in order to reduce the 
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TABLE I. Size and basic parameters of detectors along the beam. The origin of the z position corresponds to the start of 
the E391a detector. In the sizes of cross section and hole, dia. and sq. represents diameter and square, respectively. Visible 
fraction (R v i 3 ) is defined as the ratio of the energy deposit in the active material to that in the whole volume. The thickness 
is expressed in units of radiation length. The recorded number of photo-electrons per 1-MeV visible energy deposit was 10 or 
more for sandwich counters (CC02-05), and at least 0.5 for lead glass detectors (CC06-07), respectively. 



detector 


z position (cm) 


outer dimension (cm) 


inner dimension (cm) 


configuration 


Rvis 


thickness (Xo) 




zoy . i 


tjz.u Qia. 


10.0 Qia. 


lead / scint . 


U.oZ 


10. / 


CC03 


609.8 


25.0 sq. 


12.0 sq. 


tungsten/scint 


0.23 


7.6^ 


CC04 


710.3 


50.0 sq. 


12.6 sq. 


lead/scint. 


0.28 


11.8 


CC05 


874.1 


50.0 sq. 


12.6 sq. 


lead/scint. 


0.28 


11.8 


CC06 


925.6 


30.0 sq. 


15.0 sq. 


lead glass 


1.0 


6.3 


CC07 


1000.6 


30.0 sq. 


15.0 sq. 


lead glass 


1.0 


6.3 


BHCV 


1029.3 


23.0 sq. 


no 


plastic scint. 


1.0 


0.007 


BA(Run-2) 


1059.3 


24.5 sq. 


no 


quartz 


1.0 


1.5 










lead/scint. 


0.31 


13.3 


BA(Run-3) 


1059.3 


24.5 sq. 


no 


quartz 


1.0 


1.2 










PWO 


1.0 


16.8 



a The lamination of CC03 was parallel to the beam. 



7T backgrounds produced by beam-gas interactions to 
a negligible level at the sensitivity corresponding to the 
SM prediction of K\ — » ifivv. The amount of dead ma- 
terial along the path of particles from the decay vertex 
to the detector had to be minimized in order to achieve 
highly efficient detection even for low-energy particles. 
A differential pumping method was adopted for this pur- 
pose. 

The entire detection system was placed in a large vac- 
uum vessel, as shown in Fig.[5J The pressure in the outer 
part, where the detectors were placed, was around 1 Pa, 
and the pressure in the inner part, through which the 
beam passed, was 1 x 10~ 5 Pa. The two regions were 
separated by a laminated membrane sheet with a thick- 
ness of 20 mg/cm 2 . As shown in Fig.[TUl the sheet was a 
lamination of four films. The EVAL film had low trans- 
mission for oxygen gas (mostly air) and the nylon film 
strengthened the sheet. The polyethylene layers on both 
sides were used to make a tight connection by using a 
heat iron press. The bag-shaped membrane covered the 
inner surface of the Csl detectors by using a skeleton 
structure of thin aluminum pipe, similar to a camping 
tent. 

Two sets of rotary and root pump systems were con- 
nected through a manifold and eight ports to the outer- 
vacuum part. The pumping speed of each system was 
1200 m 3 /hour. Four turbo molecular pumps, each having 
the pumping speed of 800 1/sec, were connected between 
the inner vacuum part and the manifold. They produced 
the necessary pressure difference of an additional five or- 
ders of magnitude. 

For the PMTs installed in vacuum, the pressure had 
to be less than 10 Pa to prevent high voltage discharges. 
The PMT operation in vacuum additionally caused a 
cooling problem due to the absence of convection. We 
modified the configuration of the resistor chains of the 
PMTs, and cooled them with a water circulation system 
as described in Ref. [2(|. The temperature was stable 



within ±0.1° for the Csl calorimeter. 



Low density 
Polyethylene (80 |im) 




EVAL film (lSum) Nylon film (IS urn) 



FIG. 10. Membrane used for the vacuum separation. 



6. Electronics 

About 1000 PMTs were used as readout devices for all 
detectors. Signals from the detectors in the vacuum ves- 
sel were extracted through individual feed-through con- 
nectors. A connector was made of a coaxial cable simply 
molded to a metal flange with a resin. 

All detector signals were fed to Amplifier- 
Discriminator (AD) modules, which were newly 
developed for the experiment. An AD module accepted 
16 PMT signals and generated 16 analog signals, 16 
discriminated signals, and two 8-channel linear sums. 
Because early discrimination prevents time deterioration 
due to distortion of the signal shape in the cable, the AD 
modules were placed near the vacuum vessel to shorten 
the cables. 

The analog output signal was derived from each in- 
put signal with a throughput of 95%. It was sent to 
a charge-sensitive analog-to-digital converter (ADC) in 



8 



the counting hut through a 90-m coaxial cable, while the 
cable length of the other outputs was 30 m. The discrim- 
inator output was generated with a very low threshold of 
1 mV, which corresponded to an energy deposit of 1 MeV 
for the Csl calorimeter and below 1 MeV for the other 
detectors. The discriminator output was sent to a time- 
to-digital converter (TDC) in the counting hut through 
30 m of twisted-pair cable after passing through a fixed 
delay of 300 ns in the AD module. The signal summed 
over 8 inputs was sent to the counting hut through 30 m 
of coaxial cable and used to form trigger signals for data 
acquisition. 

The ADC module, LeCroy FASTBUS 1885F, had 96 
input channels, each channel having the equivalent dy- 
namic range of a 15-bit ADC in its 12-bit data by using 
a bi-lincar technique. The typical resolution at a low en- 
ergy range was 0.13 McV/bit for the Csl calorimeter, and 
0.035 McV/bit for photon veto detectors. The gate width 
for the Csl calorimeter was 200 ns. The TDC module in 
TKO (Tristan KEK Online-system),which was an elec- 
tronics platform developed in KEK [26[ , was operated at 
a full range of 200 ns and a resolution of 50 ps. The ana- 
log and discriminated signals were already delayed by 300 
ns (60-m cable with a propagation velocity of 20 cm/ns) 
compared to the timing of linear-sum signal at the en- 
trance of the counting hut. All cables to the counting 
hut were placed inside trays covered with copper-plated 
iron sheets to minimize ground noise caused by alternat- 
ing magnetic fields. The pedestal widths of almost all 
the ADC channels were less than 1 bit. 

Any failures in the PMTs lead to serious problems in 
the experiment. Almost all PMTs were operated at volt- 
ages below 60% of the rated voltage. Such a large margin 
was made possible by using very sensitive ADCs and very 
low thresholds for the TDCs. All PMTs operated prop- 
erly during the entire running time. 

Multi-hit TDCs (MTDC) were used for the BHCV and 
BA to cope with high counting rates. We did not benefit 
from using MTDCs in Run-1 because the input pulse 
width was set to 100 ns. For Run-2 and Run-3, the pulse 
width was shortened to 25 ns. 



D. Trigger and data acquisition 

The signals from eight adjacent Csl blocks (a segment) 
were collected into a summed signal; there were 72 such 
segments as shown in Fig. [TT] 

The segments were used to form a trigger signal. A 
threshold, which corresponded to an energy deposit of 
80 MeV, was applied to each summed signal, and the 
number of segments whose summed signal was above the 
threshold, Nhc, was counted. We required Nhc > 2 
for the K\ — > ifivv trigger. We also required an anti- 
coincidence of several veto counters, CV, MB, FB, and 
CC02 - CC05 with a rather high threshold. Later, in the 
offline analysis, these conditions were set tighter: a higher 
threshold for photon detection and lower thresholds for 
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FIG. 11. Segments of Csl blocks. 



the vetos. 

In addition to the trigger for — > tt°vv, we prepared 
triggers for light pulsars: xenon lamps for Csl and LEDs 
for the other detectors. They were used to monitor short- 
term drifts in the PMT gains. Moreover, they were useful 
for studying the effects of beam loading by flashing them 
within and outside the beam spills. Triggers for cosmic 
ray muons and punch-through muons were also used for 
detector calibration. 

We prepared two types of triggers to record accidental 
hits in the detector system by using an electronic pulse 
generator and the TM, which was a counter telescope 
near the target. The accidental activities observed in the 
two types of triggers were consistent with each other for 
all detectors except the BA and BHCV. While the pulsar 
trigger had no correlation with the event time and was 
randomized with respect to the asynchronous timing of 
events, the TM trigger reflected intensity variations of 
the primary beam. The TM trigger data were normally 
used because we observed the micro time-structure of the 
beam extraction. 

Accidental losses were estimated by applying the anal- 
ysis cuts used for K\ — > -K a vv to the TM trigger data. 
The estimated value was compared with the value esti- 
mated from Monte Carlo simulations with imposition of 
the TM trigger data, and consistency was confirmed. 

Data were collected through multiple parallel systems 
of the VME crates operated with a CPU to control 
them. Environmental data such as temperatures at 100 
locations, vacuum pressure at several places, and single 
counting rates of sampled channels were accumulated by 
using PCs. The CPU and PCs were distributed on the 
network. The basic software used in the E391a experi- 
ment was MIDAS j27|. 
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The dead time of the DAQ system was around 600 
/is/event. In a typical run, in which the proton intensity 
was 2.5 x 10 12 POT (protons on target) in a 2-s spill every 
4 s, the trigger rate was 300 per spill (150 Hz), and the 
live time was 91%. The data size was 3 Mbytes/spill and 
the typical data size collected per day was 60 GB. 



beam axis. Finally, the gain constants were refined by 
an iteration process based on the kinematic constraints 
of K\ — > 7r 7r°7r° decay. The short-term variation of 
the gain was corrected by using xenon-lamp light pulses. 
The reconstructed K\ mass and the width were well sta- 
bilized, as shown in Fig. IT3l 



E. Calibration 

Good linearity between the energy deposit and the 
ADC output was observed in all detector subsystems. 
Their gains were calibrated through a constant in units 
of McV/bit. The gain constants of all detectors were ba- 
sically calibrated in situ, after assembling and in vacuum, 
by using cosmic-ray muons and/or punch-through muons 
coming from the upstream region of the primary beam- 
line. While the cosmic-ray muons primarily traveled in 
the downward direction, the punch-through muons were 
parallel to the beam. In the case of sandwich detectors, 
we selected the muons so that their primary direction was 
perpendicular to the lamination. 

The Csl modules were initially calibrated by using 
cosmic-ray tracks such as the example shown in Fig. 1121 
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FIG. 12. Cosmic ray track used for calibration, 
ring shows the MB. 



The outer 



The punch-through muons were used to cross-check the 
cosmic-ray calibration because their directions were per- 
pendicular to each other and their penetration lengths 
were different (7 cm for cosmic-ray muons and 30 cm for 
punch-through muons). The gain constants were refined 
by using the two photons from a ir° produced in a spe- 
cial run in which an aluminum plate was installed on the 
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FIG. 13. Stability of the kinematic variables over the entire 
period of Run-3. The upper and lower graphs show the peak 
and width of the effective mass distribution of 3n° from the 
K% — y 7r c V°7r decay, respectively. 

In the Monte Carlo simulations, we smear ed the ene rgy 
deposit of each photon by the function a/ -J E(GcV) + b. 
The first term was given as 0.008/V^(GeV), which 
was consistent with the statistical fluctuation of pho- 
toelectron yields of 16 pe/MeV. The parameter b was 
determined to be 0.004 by tuning it to reproduce the 
K\ — > Tr Q ir°ir° invariant mass distribution. Because this 
smearing was applied to the deposited energy instead 
of the incident photon energy, the small value indicates 
good understanding of the calibration. 

The beam loading effect was examined by flashing the 
light pulse during and between the beam spills. It was 
negligibly small in the current experiment for all detec- 
tors, except for the BHCV and BA in which the PMT 
gain was shifted by 10%. 

The timings were determined relative to one of the pho- 
ton clusters in the Csl calorimeter. First, TDC constants 
(ns/count) were measured for all TDC channels by using 
a fixed delay, and the time-zero value was calibrated from 
the data. Cosmic-ray and/or punch-through muons were 
also used for the time-zero calibration. They were deter- 
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mined step-by-step by utilizing the overlapping parts of 
different detectors with respect to the muon tracks. 

The time-zero values among Csl modules were de- 
termined by using cosmic ray tracks. Finally the cal- 
ibration was refined by using six photons from Kl — > 
7r 7r 7r° decays. The time difference between two pho- 
tons was determined with a standard deviation of 0.3 ns, 
as shown in Fig. 1141 where the timing of a photon cluster 
was estimated from the timing of the central block that 
had the local maximum energy deposit. 
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FIG. 14. Time difference between two photons from the 
K L — > 7r°7r 7r° decay, obtained by subtracting the timing of 
lower energy photon from that of higher energy photon. 



III. DATA ANALYSIS 

A. Summary of three runs and outline of present 
analysis 

During the course of the experiment, there were (1 — 
2.5) x 10 12 POT per spill, with a total yield of 4.6 x 10 18 
POT for physics runs. A summary of three physics runs 
is listed in Table [XXJ The experimental setup was modi- 
fied during intervals between the runs, and the running 
parameters were refined on the basis of the results of a 
previous run. 

TABLE II. Summary of three physics runs. 



Run 


Run period 


POT 


Remarks 


Run-1 
Run-2 
Run- 3 


Feb.-Jun. 2004 
Feb.-Apr. 2005 
Oct.-Dec. 2005 


2.1 x 10 18 
1.4 x 10 18 
1.1 x 10 18 


Membrane problem 
Be absorber 
New BA 



The data quality of Run-1 was severely affected be- 
cause the membrane for vacuum separation drooped into 
the beam core, causing many beam interactions. The re- 
sults obtained from a 10% data sample of Run-1 have 



been published elsewhere [l3j . The single-event sensitiv- 
ity (S.E.S.) was (9.11 ± 0.20 stat ± 0.64 syst .) x 10~ 8 with 
1.9 ± 1.0 background events expected in the signal box 
(iVbg). The complete data sample obtained in Run-1 was 
analyzed with a blind analysis. In the blind analysis, the 
signal candidate events were not examined until after the 
final selection criteria were set, and the criteria were not 
changed later. The results of the complete data sam- 
ple of Run-1 were found to be essentially consistent with 
those obtained from the 10% data sample. The S.E.S. 
was (5.14 ± 0.25) x 10" 8 with N bg = 2.10 [H- The 
number of observed events in the signal box (N b) was 
for the 10% sample and 1 for the complete data sample. 
The corresponding upper limits for the branching ratio 
of K\ — > tt°vv mode at the 90% confidence level (C.L.) 
were found to be 2.1 x 10 -7 and 1.26 x 10 -7 , respectively. 
Because the sensitivity of the Run-1 data sample was not 
as high as those of Run-2 and Run-3, and the running 
condition of Run-1 was considerably different from the 
later runs, the results of Run-1 are not included in the 
results in this article. 

The membrane problem was rectified before Run-2. As 
described in Sec. II. C. 4, the plastic scintillator of the 
BHCV was replaced with a thicker one, and the dis- 
criminated pulse width for the BHCV and BA was re- 
duced. An additional collar counter, CC00, was installed 
in front of the FB outside the vacuum vessel. In Run-2 
and Run-3, a neutron absorber made of beryllium was 
inserted into the beam line. The results for the complete 
Run-2 data sample, S.E.S. = (2.91 ± 0.31) x 10~ 8 with 
N bg = 0.42 ± 0.14 and N ob = 0, and the upper limit of 
BR{K° L -> n°vp) = 6.7 x 10" 8 at the 90% CL, have been 
published in a letter [l4j ). 

From the previous Run-2 analysis, the dominant source 
of the background was found to be the interaction of the 
beam halo (mostly neutrons) with the detectors near the 
beam ("halo neutron background"). The halo neutron 
background was from three sources: 7r°'s produced by 
the interaction with CC02, and 7r°'s and rj's produced 
by the interaction with the CV. In the case of CC02- 
7T° background, the computed Zvtx could be shifted 
downstream due to shower leakage or photo-nuclear in- 
teractions in the Csl calorimeter. In the case of CV- 
7T° background, Zytx could be shifted upstream due to 
the fusion of multiple photons or the overlap of other 
hits in the calorimeter. In the case of CV-77 back- 
ground, Zvtx could be shifted upstream by applying the 
7T° assumption in the reconstruction of the two photons. 

In the previous Run-2 analysis, the background level 
Nbg was estimated with a different method for each back- 
ground source. A special run with an Al plate was used 
to estimate the background level for CC02-7T . A bifurca- 
tion method with a pair of cut sets was used to estimate 
the CV-7T background. An MC calculation was used to 
estimate the CV-77 background. The difference in esti- 
mating these backgrounds made it difficult to carry out 
optimizations with the objective of obtaining the best 
signal-to- noise ratio, S/N. 
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To improve the analysis of Run2-3 data as described 
in this paper, we developed a method to generate a large 
number of beam halo interactions by using the hadron- 
interaction code FLUKA [3(| , whose reliability was con- 
firmed by the data from the beam survey [I?} ■ The sim- 
ulation generated beam-halo interactions with approxi- 
mately 8 times the total number of real data obtained in 
Run-3. For the CV-77 background case, we recycled the 
event seeds to boost the r\ production and obtained ap- 
proximately 80 times larger statistics as compared to the 
Run-3 data. Although the simulation was carried out for 
Run-3, the events were used for Run-2 with minor mod- 
ifications. The running conditions were almost the same 
in Run-2 and Run-3, except that the super- layers of the 
lead scintillator sandwich of the BA were replaced with 
PWO hodoscopcs in Run-3. 

In the analysis presented in this article, we first care- 
fully checked the MC events against the data from the 
aluminum plate run to confirm the hadronic interaction 
model. Next, we optimized the selection criteria by mon- 
itoring the S /N ratio of MC data. Here, S is the number 
of events inside the signal region (340 < Zvtx < 500 cm 
and 120 < P T < 240 MeV/c) for the K° L -> tt°vD MC 
events, and N is the MC background events caused by 
beam- halo interactions [32l |. The optimization based on 
the MC is also effective to avoid potential bias in reana- 
lyzing the Run-2 data. For the data processing, we first 
masked the events in the signal region for both Run-2 
and Run-3. Once we fixed all the selection criteria by 
the optimization described above, we compared the data 
and MC for the events outside the signal region. Af- 
ter confirming the consistency, we finally examined the 
events inside the signal region. 



B. Event reconstruction 

1. Photon clustering 

The energy of each photon was measured by forming 
a cluster of Csl blocks with finite energy deposits. First, 
a search was made for a block with the locally maximum 
energy, i.e. the maximum energy among the five blocks 
geometrically sharing a side. An envelope surrounding 
the local-maximum block was developed by gathering 
blocks that shared one side and had lower energies. An 
envelope with a fixed boundary was required to reject ex- 
tra particles that hit the outsides of the two photon clus- 
ters. It was also necessary to determine that there was 
only one local maximum within the envelope, to discrim- 
inate clusters that arose from the fusion of more than one 
photon. Figure [T5l shows the distribution of the number 
of crystals contained in one photon cluster by comparing 
the six photon events of data and K\ ->■ 7r 7r°7r° MC. The 
number of crystals was obtained by counting the crystals 
having the energy deposit greater than 5 McV. The size 
of the cluster was well reproduced by the simulation. 
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FIG. 15. Distribution of the number of crystals having the en- 
ergy deposit greater than 5 MeV in a photon cluster, compar- 
ing the data and K% — > ty tt°tv MC. The top portion shows 
a comparison between the number of events in the data and 
MC, and the bottom portion shows the ratio between them. 
All analysis cuts for the Kl — > 7r 7r O 7r° mode, as will be de- 
scribed in Sec. V.B., are imposed. 



2. Energy and position correction 

The hit position of a photon at the front of the Csl 
calorimeter was first estimated by determining the cen- 
ter of energy of hit blocks in each photon cluster. The 
Zvtx coordinate was calculated by using two photon 
clusters and assuming the ir° mass. The energy, posi- 
tion, and injection angle were obtained for each photon. 

In this experiment, the energy and position of each 
photon had to be corrected for the injection angle to ob- 
tain better resolution. The thickness of the Csl crystal 
(16.2Xo) was not sufficient to prevent a small leakage of 
energy from its downstream end, and the energy leak- 
age depended on the injection angle. We accumulated a 
large sample of MC data for a single photon with various 
energies, positions, and angles, and compiled look-up ta- 
bles to correct energy and position of each photon. The 
correction was iterativcly carried out by using data from 
the tables. Figure [TrJl shows a result of the correction. 



3. Sorting of events 

The events were sorted into several samples accord- 
ing to the number of reconstructed photon clusters. The 
two-photon sample was used to search for the —> 
ifivv decay and for monitoring the —> 77 decay. 
The four-photon and six-photon samples were used for 
monitoring the K\ — > ir°ir° and — > 7r 7r°7r° decays, 
respectively. 
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FIG. 17. Schematic diagram of simulation. 



FIG. 16. Difference of the reconstructed x-position of the 
photons on the calorimeter from the true value (Ax), obtained 
by K% — > tx°uv MC simulation. The x coordinate denotes the 
horizontal position at the face of the Csl. The resolution was 
improved after the energy and position corrections. 



4- 7r reconstruction 

By using the position and energy of two photon clus- 
ters, the decay vertex (Zvtx) and the momentum of the 
7T° were reconstructed with assumptions that the invari- 
ant mass of two photons was equal to the 7T° mass and 
that the decay vertex was on the beam axis. The opening 
angle of two photons (9) was calculated from the equation 

M*„ =2-Ei£?2(l-cos0), (1) 

where E\ and E2 are the energies of the two photons. 
The momentum and transverse momentum (Pt) of the 
7r were calculated from the decay vertex. 



C. Monte Carlo simulations 

We generated Monte Carlo events that had the same 
structure as the recorded events, and analyzed them with 
the same code as for the real events. Because a full 
shower simulation would require extensive computation 
time, the calculations were separated into several stages 
and streams, as shown in Fig. 1171 The first stage was 
the beamline simulation in which the target, collima- 
tors, magnetic field, and absorbers were introduced. The 
GEANT3 package [H[ with the GFLUKA plug-in code 
for hadronic interactions was used for this simulation. 
The momentum, position in the transverse directions, 
and angle distributions at the exit of the last collima- 
tor C6 were obtained for various beam particles. 



For K^, functions for the momentum and position dis- 
tributions at the C6 exit were obtained by fitting the 
Monte Carlo data. The K\ beam beyond C6 was subse- 
quently generated with respect to the momentum spec- 
trum and targeting angle at C6. The generated K®'s 
were used for the study of various decays through 
detector simulations with the GEANT3 package. In the 
generation of the K\ — > ir^vv events, we assumed V-A 
interactions. In the detector simulations, particles were 
traced until their energy decreased below the cutoff value 
(e.g. 0.05 MeV for photons and electrons). Typical re- 
sults of Kl decay simulations: reconstructed mass, ver- 
tex position, momentum, and transverse momentum for 
the six-photon data samples are shown in Fig. 1181 Small 
corrections to the momentum and radial position of the 
K® were included. Each distribution is reproduced by 
the simulation of the if£ — > 7r 7r 7r° decays. The method 
used to reconstruct the six-photon events and the event 
selections are described in a later section. 

Another analysis stream from the MC particles at C6 
was the generation of halo neutrons. First, the core 
neutrons were removed from the neutron data sample. 
The remaining halo neutrons were used multiple times 
as seeds. During the generation, we introduced Gaus- 
sian fluctuations to the momentum, position, and angle 
of each new event to prevent duplications of the same 
event. We used the hadron code FLUKA to simulate the 
interaction between halo neutrons and the detector. Sec- 
ondary particles generated by FLUKA were fed into the 
GEANT3 detector simulation. 

The detector response was simulated by using the 
GEANT3 package. The energy deposit and the timing 
in each detector subsystem was stored. Trigger condi- 
tions were simulated according to the energy deposits of 
the detector elements, including the segments of the Csl 
shown in Fig. [TT] 
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FIG. 18. Distributions of reconstructed mass (top left), vertex position (top right), momentum (bottom left), and transverse 
momentum (bottom right) obtained from the six-photon data samples, compared to the data and the K° L — ► 7r°7r°7r MC. In 
each plot, the top portion shows a comparison between the number of events in the data and MC, and the bottom portion 
shows the ratio between them. All the veto cuts and kinematic selections, except those for the respective abscissa variable for 
the K° L — >• 7r c V 7r° mode, are applied in these plots. Details of the cuts are described in Sec. V.B. 



D. Reproducibility of the MC simulations 

In addition to the distributions of kinematic variables 
for the Kl —> 7r 7r°7r° decay mode, the events with four 
photons reconstructed in the calorimeter were analyzed 
in order to verify the detection inefficiency of photon 
counters in the simulation. Figure Q1J] shows the recon- 
structed mass distribution for four-cluster events. In ad- 
dition to the Kl —> 7r°7r° events at the K\ mass, there 
was a tail in the lower mass region due to contamination 
from Kl — > 7r^7r^7r^, two out of six photons of which es- 
caped detection. The number of events in the tail was 



reproduced by our simulations. 

A simulation with the FLUKA hadronic-intcraction 
package to a dedicated run ( "Al plate run" ) was carried 
out in order to confirm the reproducibility of hadronic in- 
teractions. In the Al plate run, a 0.5-cm-thick aluminum 
plate was inserted into the beam at 6.5 cm downstream 
of the rear end of CC02, as shown in Fig. [501 Because 
the position of the Al plate is known, the invariant mass 
of two photons can be reconstructed from the energy and 
position of two photons. Figure [5T] demonstrates the sim- 
ulation in which the invariant mass distribution (from 
7T° mass to rj mass) of the events from the Al plate run 
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FIG. 19. Reconstructed invariant-mass distribution of events 
with four photons in the calorimeter. The points show the 
data and the histograms indicate the contribution of if£ — >■ 
7r°7r° and iv£ — > tt°ty°tt decays (and their sum), as expected 
from the simulation, normalized by the number of events in 
the Kl — > tt°ty° peak. All veto and kinematic cuts except for 
the cut on the 47 invariant mass for the K% — > 7r°7r° mode are 
applied. It includes the shower-leakage correction, but does 
not include the correction for photo-nuclear interactions. 



was reproduced with two photons in the calorimeter. 

Next, by assuming the ir° mass and reconstructing 
the vertex position in the Al plate run, the effect of 
shower leakage and photo-nuclear interactions in the 
Csl was examined. Figure [22] shows the distribution of 
Zvtx obtained from the Al plate run by comparing the 
data and MC. If the photon energy is incorrectly mea- 
sured as being low due to these processes, the ir° vertex 
position is reconstructed downstream of the true position. 
Because the true ^vtx is known in the Al plate run, we 
can estimate the effect of these processes by using the 
tail in the distribution of the reconstructed Zvtx- 

Photo-nuclear interactions arc not implemented in 
GEANT3, and therefore the effects were estimated from 
GEANT4 simulations [H]. We estimated the probabil- 
ity of photo-nuclear interactions by taking the difference 
between MC samples with and without photo-nuclear in- 
teractions. The energies of two photons were smeared 
according to the probability of the interaction. As shown 
in Fig. [22j a tail on the downstream side of the peak is 
reproduced by the simulations with photo-nuclear inter- 
actions. 



IV. EVENT SELECTION AND BACKGROUND 
ESTIMATION 

A. Event selection 

Event selection consisted of veto cuts and kinematic se- 
lections. Vetoing extra activities was the primary method 
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FIG. 20. Schematic layout of the Al plate run. 
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FIG. 21. Reconstructed invariant mass distribution of the 
two-photon event in the Al plate run with all the veto cuts 
imposed. The points represent data and the solid line is from 
FLUKA simulations. The peaks at 0.14 GeV/c 2 0.55 GeV/c 2 
correspond to 7r°and r\ particles, respectively. 
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FIG. 22. Reconstructed Zvtx distribution of two photon 
events in the Al plate run. The points represent the real 
data and the solid (dashed) line is obtained from the FLUKA 
MC simulaiton with (without) the effects of photo-nuclear 
interaction. In both cases, the shower-leakage correction was 
applied. 
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used to isolate the — > tt°vv signal events from the 
background events. Kinematic selections were applied to 
obtain further rejection of background events. 



1. Veto cuts 

Particle veto with a hermetic detector system is the 
primary method to reject possible background sources. 
For example, events made by the if£ — > 7r°7r° or — >• 
ir + ir~ir° mode can be rejected by veto cuts because these 
modes have at least two extra activities. The method also 
works for halo neutron backgrounds because the hadronic 
interactions of halo neutrons are often accompanied by 
extra particles such as protons or pions. 

Energy thresholds for veto cuts were set at two lev- 
els. For loose veto thresholds (typically, 10 MeV), no 
timing cut was made to reduce inefficiencies caused by 
accidental hits before the event time. For tighter veto 
thresholds (typically, 1-2 MeV), events with a TDC value 
within roughly ±5er from the on-time peak were rejected. 
Combining these conditions, we balanced the efficiencies 
of detecting extra particles and the suppression of the 
acceptance loss caused by accidental hits. 

Csl veto cut: In addition to its main role of mea- 
suring the energy and position of photons, the Csl also 
served as a veto detector for extra photons. Extra activ- 
ities that were not reconstructed as photon clusters, i.e. 
"single crystal hits," were rejected with this cut. How- 
ever, a photon occasionally creates a single crystal hit 
near its genuine cluster due to fluctuations in the elec- 
tromagnetic processes. Thus, applying a tight energy 
threshold for a single crystal hit near the photon cluster 
can cause acceptance loss for signal candidates. To re- 
cover this loss, the energy threshold for a single crystal 
hit was determined as a function of the distance d to the 
closest cluster as shown in Fig. [531 



• Sthrcs. = 10 MeV for d < 17 cm, 

5 - (3/8) (d - 17) MeV for 17 < d < 25 



• -Ethre; 

cm, 



E 



thr 



2 MeV for d > 25 cm, 



where -Ethres. is the energy threshold for a single crystal 
hit. 

Other photon veto cuts: Photon veto detectors con- 
sisted of barrel counters (MB and FB), beam collar coun- 
ters (CC00, CC02 - CC07), and a beam hole counter 
(BA). 

Because the MB was read out at both the upstream 
and downstream ends, the energy deposit was determined 
as the geometrical mean of the visible energies in both 
sides in order to cancel the position dependence of the 
light yield. The timing requirement for the MB was loos- 
ened when a tighter energy threshold was used because 
the timing was broadly distributed in the halo neutron 
backgrounds. When a photon hit the Csl calorimeter, 
low-energy photons and electrons that were created in 
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FIG. 23. Energy deposition in a single crystal hit versus the 
distance from the nearest photon cluster of the two-photon 
events from the K% — > ty°vu (top) and K% — ► 7r°7r° (bottom) 
simulations. Events above the solid line were rejected. 



an electromagnetic shower occasionally went backward 
and hit the MB. This process caused a larger acceptance 
loss in the MB as compared to the loss in other veto 
detectors. 

For the collar counters in the downstream region 
(CC06 and CC07) , both tighter and looser energy thresh- 
olds were set higher than for other detectors, in order to 
reduce accidental losses due to beam-induced activities. 
The BA located inside the beamline rejected photons es- 
caping into the beam hole. The energy thresholds for 
these counters are summarized in Table IIIII 

Charged particle veto cuts: Background events 
produced by charged decay modes were mostly re- 
jected by cuts on the charged particle veto detectors, 
mainly CV, BCV, and BHCV. As in the case of MB, 
an energy deposit in the BCV was determined as the 
geometrical mean of the hits in the upstream and down- 
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stream ends, and the timing requirement was loosened. 
The energy thresholds for these counters are summarized 
in Table [ 



TABLE III. Summary of tighter veto energy threshold for 
each detector. 



Detector 


Threshold 




1 IMC V 


r J_> 


1 A/T^V a 

i ivie v 


ppno 


i ivie v 


iiOV 


U. ( o Me V 


MB Inner 


1 MeV b 


MB Outer 


1 MeV b 


CV Outer 


0.3 MeV 


CV Inner 


0.7 MeV 


CC03 


2 MeV 


Csl 


d6p6nds on d c 


Sandwich 


2 MeV 


CC04 Scintillator 


0.7 MeV 


CC04 Calorimeter 


2 MeV 


CC05 Scintillator 


0.7 MeV 


CC05 Calorimeter 


3 MeV 


CC06 


10 MeV 


CC07 


10 MeV 


BHCV 


0.1 MeV 


BA Scintillator (Run-2) 


20 MeV^ 


BA PWO (Run-3) 


50 MeV 3 


BA Quartz 


0.5 MIPs"3 



a Sum of inner and outer layers. 

b Geometrical mean of upstream and downstream ends. 
c Detailed in the text. 
d Summed over layers. 

e Determined by AND logic of scintillator/PWO and quartz. 



2. Kinematic selections 

Kinematic selections were applied to discriminate the 
signal from background events by using the information 
of two photons in the Csl calorimeter. These selections 
can be categorized into three types: photon-cluster qual- 
ity selections, selection on photons, and ir° selections. 

Photon-cluster quality selections: Photon-cluster 
quality selections were developed from the information 
of each photon cluster. The energy was required to be 
E H > 250 MeV and E L > 150 MeV, where E H and E L 
are the energies of the higher- and lower-energy photons, 
respectively. We also made selections on the incident po- 
sition, size, profile of the energy distribution, and timing 
dispersion of each photon cluster. These cuts were ap- 
plied to ensure that two photon clusters were from the 
A'£ — } 7T° vv decay and were distinguished from fake clus- 
ters created by hadronic showers or photons produced by 
halo neutron interactions. 

In addition, two neural network selections were devel- 
oped pertaining to the shape of a photon cluster. The 
first neural network was used to reduce events with clus- 



ters that overlapped those from other photons or associ- 
ated particles and reconstructed as a single cluster ( "fu- 
sion cluster"). In the case of K\ — > ir°ir° background, a 
fusion cluster results in an inefficiency to extra photons; 
in the case of CV-7T background, it results in a larger 
photon energy such that the event vertex comes into the 
signal region. This neural network was trained by single 
and fusion clusters obtained from the Ay — > A MC 
samples. The second neural network selection was used 
to reject the CV-77 background. Because 77 particles were 
produced near the front end of Csl and had a larger in- 
variant mass than a 7r°, the two photons produced by 
77 decay tended to have large incident angles to the Csl 
calorimeter. This neural network was trained by the MC 
samples of A£ ir°vv and CV-77 background. The re- 
jection power of these two neural network selections is 
shown in Fig. 1241 with the distribution of the output 
values from the neural network for the signal and back- 
ground modes. 
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FIG. 24. Distributions of the output values obtained by two 
neural network selections for the signal mode (hollow) and 
the background events (hatched). The top figure shows the 
fusion neural network selection, in which the background is 
CV-7T . The bottom figure shows the CV-77 neural network 
selection, in which the background is CV-77. In both plots, all 
analysis cuts except for the respective neural network cut are 
imposed for the signal mode. Several cuts are loosened for 
the background case in order to enhance the events. 

Selection on photons: We required the distance be- 
tween two photons to be greater than 15 cm, and the 
timing difference between them to be —9.6 < T(Eh) — 
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T(E L ) < 18.4 ns, where T(E H ) and T{E L ) are the tim- 
ing of the higher and lower energy photons, respectively. 
In addition, the energy balance of two photons defined 
as {Eh — El) /{Eh + El) , was required to be less than 
0.75. These selections were needed mainly to suppress 
accidental activities. 

7T° selections: By using the information for recon- 
structed 7r°'s, several selection criteria were imposed. We 
required the kinetic energy of the reconstructed 7r°'s to 
be less than 2 GeV to suppress neutron backgrounds with 
high energy. In addition, the reconstructed 7r°'s had to be 
kinematically consistent with a K® — > n up decay within 
the proper K® momentum range. A neural-network se- 
lection was applied to the discrepancy between two an- 
gles. One angle was calculated by connecting the incident 
position of photons and the position of reconstructed 
7r°'s, and the other angle was estimated by the neural 
network based on the cluster shape. This selection sup- 
pressed CV-7T and CV-ij backgrounds because the dis- 
crepancy of the angle increased in these processes. To re- 
duce the Kl — >• 77 background, we calculated the open- 
ing angle between the two photon directions projected 
onto the Csl calorimeter plane, and required it to be less 
than 135 degrees ( "acoplanarity angle cut"). 



3. Signal region 

We set the signal region in the Zvtx-Pt plane to be 
340 < Z V tx < 500 cm and 0.12 < P T < 0.24 GcV/c. 
The requirement on ^vtx was determined to avoid back- 
ground 7r°'s coming from the CC02 (upstream) and CV 
(downstream). The lower limit for Pt was set to reduce 
the K\ — > 7r°7r° and CV-ij backgrounds; these events 
cluster in the low Pr region due to kinematics and a 
large halo neutron flux near the beam center, respec- 
tively. The upper bound for Pt was determined by the 
kinematic limit of the — > tt°pp decay, whose maxi- 
mum value is 0.231 GeV/c. 



B. Background estimation 

1. Halo neutron background 

The halo neutron background was estimated from com- 
bined MC simulations with FLUKA and GEANT3. 

Figure [231 shows the distribution of background events 
in the Zvtx-Pt plane as estimated from the MC sam- 
ples. The outside of the signal box was divided into four 
regions, as shown in Fig. 1251 The number of events in 
each region was dominated by the halo neutron back- 
grounds. The numbers of halo neutron backgrounds were 
normalized to the number of events in the upstream re- 
gion (shown as Region- (1)). The numbers for events out- 
side the signal box were compared between the data and 
the MC, as listed in Table IIV1 The numbers of events 



in each region are consistent within the statistical uncer- 
tainties between the data and the MC. Possible impacts 
of discrepancies outside the signal box on the estimation 
inside the signal box are included in the systematic un- 
certainties below, and will be discussed in detail in a later 
section (Sec. VLB.). 

For estimating the CC02-7T background, the effects of 
shower leakage and photo-nuclear interactions described 
in Sec. III.D were taken into account. After applying all 
event selections, we estimated 0.66 ± 0.33 sta t ± 0.20 syst 
background events in the signal box for the combined 
data sample of Run-2 and Run-3. 

For CV-7T events, a fusion cluster or a cluster made by 
a neutron produces background. After applying all event 
selections, no events inside the signal box were obtained 
with 3 times larger statistics as compared to the com- 
bined data sample of Run-2 and Run-3. An upper limit 
of 0.36 events at the la level was set for this process. 

CV-77 events were strongly suppressed by the dedicated 
neural- network cut on the cluster shape. A value of 0.19± 
0-08 s tat ± 0.10 S y S t events was estimated for this process. 
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FIG. 25. Scatter plot of Pt vs. the reconstructed Z position 
after applying all cuts to the MC samples for background sim- 
ulation, including halo neutron and K% — > 7r°7r° backgrounds. 
The region bounded by the solid line shows the signal box. 
Events around Zvtx = 275 cm are halo neutron events 
reconstructed at the position of CC02, and those around 
Zvtx = 560 cm are events reconstructed at the CV. The 
numbers of events in each region are listed in Table IIVI 



2. K L background 

The K\ background was estimated by using GEANT3- 
based MC simulations that were performed separately for 
Run-2 and Run-3 conditions. Among the backgrounds 
from Kl decay, the K\ — » 7r 7r mode made the largest 
contribution. The amount of Monte Carlo statistics for 
the K\ — >• 7r°7r° mode was roughly 70 (60) times that for 
the real data of the Run-2 (Run-3) sample. In Run-2, 
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TABLE IV. Estimated numbers of events outside the signal 
region, compared with the combined data samples of Run-2 
and Run-3. Note that the numbers of events estimated by the 
MC simulation were normalized in Region- (1). Errors in the 
MC estimates indicate statistical uncertainties. 



TABLE V. Estimated number of background events. 



Region 



Data MC estimation 



Region-(l) (CC02) 360 360.0 ± 15.6 

Region-(2) (CV) 101 77.2 ± 5.6 

Region-(3) (upstream) 8 5.9 ±1.1 

Region-(4) (low-P T ) 8 2.9 ± 0.9 

Signal box (masked) 0.87 ± 0.34 



Background source 


Estimated number of BG 


halo neutron BG CC02-7r u 


0.66 ±0.39 


CV-7T 


< 0.36 


CV-7? 


0.19 ± 0.13 


Kl decay BG -> 7r u 7r u 


(2.4 ± 1.8) x 10 -2 


K\ -> 77 


negligible 


charged modes 


negligible (0(nr 4 )) 


other BG backward n" 


< 0.05 


residual gas 


negligible (G>(HT 4 )) 


total 


0.87 ±0.41 



two events remained after applying all event selections in 
the simulations, and in Run-3, no events remained. The 
two remaining events both corresponded to the case that 
two photons in the Csl came from the same 7r° . One 
extra photon with high energy (~ 1 GeV) went through 
the Csl ("punch-through"), and the other extra photon 
with low energy (~ 10 MeV) hit the MB and failed to be 
detected by the fluctuation in an electromagnetic shower. 
The estimated number of background events from — » 
7r°7r° decay was estimated as (2.4±1.8 8ta t±0.2 sys t) x 10~ 2 
for the combined data sample of Run-2 and Run-3. 

For the MC if£ — > 77 samples whose statistics cor- 
responded to 3.8 and 8.7 times the Run-2 and Run- 
3 data, respectively, all of the event selection cuts ex- 
cept for the acoplanarity angle cut were imposed and no 
events remained in the signal region. Furthermore, the 
acoplanarity angle cut strongly suppressed the — > 
77 events, with a typical rejection of 1 x 10 5 . We con- 
cluded that the background due to K\ — > 77 was negli- 
gible in this analysis. 

For charged modes, the branching ratios of the K® — > 
7T + 7r~7r° and K\ — > irlv (I — e,/i) modes are too 
large to generate sufficient statistics in the simulations. 
The background from these modes were estimated by 
assuming that the inefficiency of charged-particle ve- 
tos was 1.0 x 10" 4 for the CV, and 1.0 x 10" 3 for 
BCV and BHCV 0. By applying the event weight 
calculated from the inefficiency, the numbers of back- 
ground events from these modes were estimated to be 
4.2 x 10" 4 for K° L -> ^e+v and less than 1.0 x 10~ 4 
for K\ — > 7r + 7r~7T°. Here, the — ► ir~e + i/ mode had 
the largest contribution to the background events among 
K\ — > irlv modes because a charge-exchange interaction 
of a 7r~ (ir~ +p — > 7T° + n) and annihilation of e + prevent 
the detection of charged particles. 



stream direction. Because we were unable to discrimi- 
nate whether photons came from the front or the back, 
these events were reconstructed as coming from the front 
and ended up inside the signal box. We estimated this 
background by using FLUKA simulations with 20 times 
the amount of statistics as compared to the combined 
data sample of Run-2 and Run-3. No events remained 
after applying all cuts, and the number of background 
events was determined to be less than 0.05 events. 

The second additional background source is the resid- 
ual gas background that occurs from interactions of beam 
neutrons with the residual gas. This process is well sup- 
pressed by a high-vacuum system that provided a very 
low pressure of 10 -5 Pa. To estimate the background, 
we carried out a dedicated run at atmospheric pressure, 
accumulating statistics that were roughly equivalent to 
0.6% of the combined data sample of Run-2 and Run-3. 
We obtained 6867 candidate events with loose event se- 
lections. By considering a reduction factor of 10 -10 due 
to the air pressure, this background was concluded to be 
negligible. 

We investigated possible contributions from unknown 
background sources by loosening the selection criteria for 
the data and MC. The numbers of events in the data sam- 
ple with loosened cuts were consistent with the prediction 
with the known background sources. 

In total, the estimated number of background events 
was 0.87±0.41, where the CV-7T and the backward-going 
7T° backgrounds were not included. Table [V] summarizes 
estimates for the numbers of background events. 



V. SENSITIVITY AND RESULTS 
A. Principle of normalization 



3. Other background sources 

In addition to the backgrounds from halo neutron 
and K\ decay, two other possible background sources 
were considered. The first one is the backward-going 
7T° background. When the halo neutron interacts in the 
end cap of the vacuum vessel located 2 m downstream 
of the Csl, 7r°'s are occasionally produced in the up- 



The single event sensitivity (S.E.S.) is represented as 

1 



S.E.S.(K° L -> = 



N(Kl decays) x A signal 



(2) 



where N(K® decays) denotes the number of K^'s that 
decayed in the fiducial region and ^4 s i gna i denotes the ac- 
ceptance for the signal mode. The value of AT(if£ decays) 
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is determined from the normalization mode as 

jxrdata 

N{Kl decays) = f^- , (3) 

where A^orm is the number of observed events from the 
normalization mode and A noTm . and BR noTm . represent 
the acceptance and branching ratio of the mode, respec- 
tively. Substituting Eq. [3] into Eq. dj the S.E.S. is ob- 
tained from the number of remaining events in the nor- 
malization mode and the ratio of acceptances between 
the signal and the normalization modes. By taking this 
ratio, uncertainties arising from variations in beam con- 
dition, etc., can be canceled. 

We examined three decay modes, namely, K® 
7r°7r°7r , K° L -> 7r°7r°, and K\ -> 77 , to obtain the num- 
ber of K\ decays, because they were fully reconstructed 
and clearly identified. Because these three decay modes 
have different numbers of photons in the final states, we 
could cross-check the reliability of the clustering method. 
The mean value of the accepted K® momentum was also 
different, which provided a check in the full range of the 
accepted K\ momentum for the K\ — > tt°vv mode. 



For the A"£ — > 77 mode, if£'s were reconstructed from 
the two photons in the Csl calorimeter by assuming the 
K\ mass. 

Several analysis cuts were imposed on the events in 
each decay mode to remove contaminations of other de- 
cay modes. Photon veto cuts were important to detect 
the modes with extra photons, such as K® L — > 7r 7r°7r° to 
the K° L -> 7T 7r° events, and K° L ->■ 7r°7r° to the 
K\ — > 77 events. The energy threshold for the veto 
in each subsystem was the same as that used in the 
K\ — > it I'D analysis (see Table IPTI]) . The reconstructed 
mass of six photons in the K a L — > 7r 7r°7r° events and four 
photons in the — > ir°ir° events had to be consistent 
with the K° L mass (from 0.481 to 0.513 GeV/c 2 ). The 
decay vertex point of K\ also had to be located in the 
fiducial region (from 340 to 500 cm), as in the case of 
K\ — > ifivv. Some additional cuts, such as the neural 
network fusion cluster selection, were also imposed on 
these decay modes. 



2. Number of Kl decays 



B. Analysis of normalization modes 



For the K° 



1. K% reconstruction 



7r 7r°7r and K° 



were reconstructed from the six and four photons in the 
Csl calorimeter, respectively. In the reconstruction, the 
number of possibile combinations of pairs of the two pho- 
tons was ( 6 C* 4 X4C2 x 2 C 2 )/3! = 15 for K° L -> 7r 7r°7r° and 
(4C2 x 2 C* 2 )/2! = 3 for K° L -> 7r°7r°. The decay vertex lo- 
cation can be calculated from the energy and position of 
the two photons of each pair by assuming the ir° mass. 
To find the best combination of the photons, the vari- 
ance in the reconstructed vertex points was calculated, 
named "pairing Xz" f° r an the possible combinations. 
The pairing \z was defined as 



E 



(Zt - zf 



(4) 



where i runs over the two-photon pairs reconstructing 
tt°, (e.g., n = 3 for K° L 7r 7r°7r° and n = 2 for K° L -> 
7r 7r .), Zi is the vertex point of i-th two-photon pair, Ui 
is the resolution in reconstructing Zi calculated from the 
energy and position resolutions of the two photons, and 



(5) 



The decay vertex of the K\ was determined as Z for 
the combination with minimum \ z . We required the 
minimum pairing \\ ^° ^ c ^ css than 3.0 and the difference 
between the next-to-minimum one to be greater than 4.0, 
in order to reduce incorrect paring. 



The acceptance of each mode was estimated from 
Monte-Carlo simulations. Because, in the simulations, 
K^s were generated at the exit of the last collimator 
(C6), the probability of decay in the fiducial region, 
(2.14 ± 0.02)%, was calculated separately and was taken 
into account in calculating the number of K® L decays. 
Losses due to accidental activities in the detector were 
also included in the acceptance. 

Table I VII shows a summary of the estimated numbers 
of K\ decays from the three decay modes K\ — > 7r 7r°7r°, 
K\ — > 7r°7r°, and K\ — > 77. The difference among the 
three modes was within the systematic uncertainties, and 
considered to come from the Csl veto, because it has the 
largest systematic uncertainties as well as a dependence 
on the number of photons. We adopted the number ob- 
tained from K° L 7T 7r°, (8.70 ± 0.17 sta t ± 0.59 syst ) x 10 9 
for the combined data sample of Run- 2 and Run-3, as the 
normalization of this analysis because the energy distri- 
bution of photons in the Csl calorimeter from K Q L — > 
tt°tt° was similar to that expected to K® L — > ifivv. Esti- 
mates of the systematic uncertainties are described later. 



C. Acceptance and single event sensitivity 

1. Signal acceptance 

The acceptance for the — > tt°vv decay was esti- 
mated from Monte-Carlo simulations. The raw accep- 
tance was calculated by dividing the number of remaining 
events after imposing all the analysis cuts by the number 
of — > tt°vv decays generated in the simulation. It 
was 1.40% for the Run-2 and 1.39% for the Run-3 data, 
respectively 
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TABLE VI. Estimated numbers of Ki decays calculated from the three decay modes in the combined sample of Run-2 and 
Run-3 data. Uncertainties in the acceptances are statistical ones due to the amount of the MC samples. N^im is the number 
of events obtained from the real data after imposing all the analysis cuts. In the _ff£ — >• 7r°7r° mode, N^otm IS obtained by 
subtracting the contamination from the K% — ► 7r°7r 7r° mode. Statistical uncertainties in the number of K% decays include 
an ambiguity in iV^, and systematic ones include an ambiguity from the reproducibility of the MC (described later) and 
statistical uncertainties of the MC in the acceptance estimate. 

mode acceptance ~Nnmm A f (_R"£ decays) 

Kl ->• 7r u 7r u 7r u (7.21 ±0.06) x lO -5 118334 (8.41 ± 0.03 sta t ± 0.53 syst ) x 10 a 
K% -> 7r°7r° (3.42 ± 0.03) x 10~ 4 2573.9 (8.70 ± 0.17 sta t ± 0.59 syst ) x 10 9 
Kj -»■ 77 (7.18 ± 0.03) x 10~ 3 35367 (9.02 ± 0.05 st at ± 0.51 ByBt ) x 10 9 



In the raw acceptance, the losses caused by geometrical 
acceptance, veto cuts, kinematic selections, and selection 
on fiducial region were included. The geometrical accep- 
tance to detect two photons in the Csl was approximately 
20%. The loss by veto cuts, or self-vetoing, was about 
50%, which was dominated by the Csl and MB vetoes. 
The loss in the Csl veto was caused by hits accompa- 
nied by the genuine photon cluster, as shown in Fig. 1231 
and the loss in the MB was caused by photons escap- 
ing from the front face of the Csl and hitting the MB. 
The acceptance of the kinematic and signal region se- 
lections was approximately 15%. These values were esti- 
mated through MC simulations of the K® L — > ir°vv mode. 
The evaluation of the acceptance loss was supported by 
the fact that the acceptance losses in the normalization 
modes {K° L -» 7r°7r 7r°, K° L -> tt°tt°, and K° L -> 77) were 
reproduced by the simulation. 

The acceptance loss due to accidental activities in the 
detector was estimated from real data taken with the TM 
trigger. The accidental loss was estimated to be 20.6% 
for the Run-3 data, in which the losses in MB (7.4%) and 
BA (6.4%) were major contributions. For Run-2, the ac- 
cidental loss was estimated to be 17.4%; the difference 
between Run-2 and Run-3 was due to the difference in 
the BA counters used in the data taking. The acceptance 
loss caused by the selections on the timing dispersion of 
each photon cluster and on the timing difference between 
two photons was estimated separately by using real data 
and was obtained to be 8.9%. Thus, the total accep- 
tance for the K° L -> ix^vv was (1.06 ± 0.08)% for Run-2 
and (1.00 ± 0.06)% for Run-3 case, where the errors are 
dominated by the systematic uncertainties that are dis- 
cussed later. Figure [26] shows the distribution of the MC 
K\ — > ir^vv events in the scatter plot of Pt-^vtx after 
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FIG. 26. Density plot of Pt vs. the reconstructed Z position 
for the K% — ¥ tPvv Monte Carlo events after imposing all 
of the analysis cuts. The box indicates the signal region for 



D. Results 



After finalizing all of the event selection cuts, the can- 
didate events inside the signal region were examined. No 
events were observed in the signal region, as shown in 
Fig. 1271 An upper limit for the — > iftvv branching 
ratio was set to be 2.6 x 10 -8 at the 90% confidence level, 
based on Poisson statistics. The result improves the limit 
previously published [lH by a factor of 2.6. 



imposing all of the other cuts. 



2. Single event sensitivity 

By using the number of K\ decays and the to- 
tal acceptance, the single event sensitivity for K Q L — ¥ 
n°uD was (1.84 ± 0.05 sta t ± 0.19 syst ) x 10~ 8 for Run- 
2, (2.80 ± 0.09 sta t ± 0.23 syst ) x 10~ 8 for Run-3, and 
(1.11 ± 0.02 stat ± O.lOsyst) x 10" 8 in total. 



VI. SYSTEMATIC UNCERTAINTIES 

Although the systematic uncertainties were not taken 
into account in setting the current upper limit on the 
branching ratio, we will describe our treatment for them 
to provide a thorough understanding of the experiment. 
In particular, systematic uncertainties of the single event 
sensitivity and background estimates due to halo neu- 
trons arc discussed in order. 
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FIG. 27. Scatter plot of Pt vs. the reconstructed Z position 
for the events with all of the selection cuts imposed. The box 
indicates the signal region for A'£ — > ifivv. 



A. Uncertainty of the single event sensitivity 

The systematic uncertainty of the single event sensi- 
tivity was evaluated by summing the uncertainties of 
the number of K\ decays and the acceptance of the 
K a L — > tt q i/D decay. Because the calculation of the former 
also includes the acceptance of the normalization modes, 
the acceptance of both the normalization mode and the 
K\ — > ifivv mode are relevant to the acceptance evalu- 
ation by the Monte Carlo simulations. To estimate the 
uncertainties in the acceptance calculation, we utilized 
the fractional difference between data and the simulation 
in each selection criterion, defined by the equation 



F' 



A 



data 



-4 



MC 



.4 



(6) 



data 



where A data and A MC 



denote the acceptance values of 
the «-th cut, calculated as the ratio of numbers of events 
with and without the cut, for the data and MC simula- 
tions, respectively. In F 1 , the acceptance was calculated 
with all the other cuts imposed. The systematic uncer- 
tainty of the acceptance was evaluated by summing all 
the fractional differences in quadrature, weighted by the 
effectiveness of each cut, as 



X/i=cuts(-^ 1 A^data) 



syst. 



-cuts (l/4ata) 2 



three decay modes used in the 



(7) 



7r°7r°, 



and K° L 



normaliza- 
— > 77, the 



For the 
tion, K\ 

calculated uncertainties were 5.2%, 5.7%, and 3.6% (in 
Run-3), respectively. The acceptances of the Csl veto 
cut had the largest uncertainties in all decay modes. 
The number of K° decays was obtained by using the 
A'£ — > 7r°7r° mode, and thus its uncertainty is quoted as 
the same value for the K\ — > ir°ir° mode. 



For the acceptance of the — > n iyp mode, the same 
systematic uncertainties as the K\ — > ir°ir° mode were 
adopted because there were no signal candidates in the 
data to be compared with the MC simulations. 

The systematic error of the single event sensitivity was 
evaluated to be a quadratic sum of the uncertainties of 
the number of K\ decays and the acceptance of the 
K° L ->• tt°i>D decay. It was 10.3% in Run-2 and 8.2% 
in Run-3, respectively. 



B. Uncertainty of the halo neutron backgrounds 

The systematic errors of halo neutron backgrounds 
were also estimated by utilizing fractional differences be- 
tween data and the simulations. There were two kinds 
of cuts in our analysis, veto cuts and kinematic selec- 
tions. For the former, the method was unable to be used 
directly because veto cuts had been applied in the early 
stages of the MC simulations to save the computing time. 
Thus, for the systematic uncertainty due to veto cuts, the 
same value obtained in the — > tt°tt° analysis was as- 
signed. For the kinematic selections, the same method 
as described in the previous section was used. The ac- 
ceptance of each cut was calculated with all veto cuts 
imposed except kinematic selections. In total, the sys- 
tematic uncertainties due to fractional differences were 
calculated to be 31%, 32%, and 44% for CC02-7T , CV- 
7T°, and CV-77 backgrounds, respectively. 

In addition, the uncertainties in the normalizations of 
the MC simulations were taken into account. Because 
the normalization was determined by using the num- 
ber of events in the CC02 region (Region-(l)), there can 
be ambiguity in estimating the CV-related backgrounds. 
As shown in Table IIV1 there was a 24% difference be- 
tween data and MC simulation in the downstream re- 
gion (Rcgion-(2)). Even though they were statistically 
consistent, we assigned the difference as an additional 
systematic uncertainty of CV-7T . For the CV-r; case, 
further ambiguity due to the reproducibility of r\ produc- 
tion should be considered. It was estimated to be 24% 
from the difference between data and MC simulations in 
the numbers of 77 events in the Al plate run, as shown in 
Fig.dU 

The total systematic uncertainties were calculated by 
summing up the contributions quadratically, to be 31%, 
40%, and 55% for CC02-7T , CV-tt , and CV-r) back- 
grounds, respectively. 



VII. CONCLUSION AND DISCUSSION 

The E391a experiment at the KEK-PS was the first 
dedicated experiment for the — »• tt°vD decay. Com- 
bining the periods of Run-2 and Run-3, the single event 
sensitivity reached 1.11 x 10 -8 . No events were observed 
inside the signal region and the new upper limit on the 
branching ratio of the K\ — > ifivv decay was set to be 
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BR < 2.6 x 1CT 8 at the 90% confidence level. The result 
improves the previous published limit given by the Run-2 
analysis by a factor of 2.6, and the E391a experiment as a 
whole has improved the limit from previous experiments 
by a factor of 20. 

The E391a experiment was also the first step of our 
step-by-stcp approach toward the accurate measurement 
of the K® — > tjPvv decay. The experiment was planned 
to confirm our experimental approach, and this purpose 
was well achieved. Several points need to be noted. 

First, we found solutions to several technical issues, 
such as the pencil beamline, differential pumping for 
ultra-high vacuum, low-threshold particle detection with 
a hermetic configuration, in situ calibration, etc., which 
can be successively used in the next step. We encoun- 
tered some technical problems that exist in the current 
apparatus, such as insufficient thickness and segmenta- 
tion of the Csl calorimeter, the structure of the CV, 
the limitation of the BA in an environment with higher 
counting rate, etc.. They will be improved in the next 
experiment. 

Secondly, we were able to control the systematic un- 
certainties to be small in the estimate of the single event 
sensitivity, as described in the previous section. A small 
systematic error is essential for an accurate determina- 
tion of the branching ratio. We are confident that the 
branching ratio of the decay — > ifivv can be mea- 
sured accurately with this method. 

A third point concerns the understanding and estima- 
tion of backgrounds. In the experiment for the decay 
K\ — > 7T° vv, elimination of all possible backgrounds is 
the only effective way to identify the decay, and it can be 
achieved by a profound understanding of backgrounds. 

The dominant source of backgrounds was the result 
of beam interactions. Although the background from 
this source was considered to be as serious as those 
from A° decays in the experiment, its understand- 
ing and estimation were very difficult to assess before 
the work reported here. The background mechanisms 
were clearly understood and divided into three different 
sources, CC02, CV-tt , and CV-ry. Methods were de- 
veloped to estimate them with rather small systematic 
errors. Based on this experience, we clearly know the 
direction of upgrades to minimize the beam backgrounds 
in the next step. 

The backgrounds from other A£ decays were reduced 
to be negligibly small in the current experiment, mainly 
due to the success of applying tight vetoes with a her- 
metic configuration. One of the important results is the 
invariant mass distribution of the events with four pho- 
tons, shown in Fig. 1191 The distribution was well repro- 
duced by the simulations. In particular, the low-mass 
region was well described by A£ — > tt tt tt decays with 



a small contamination by mis-combination events of four 
photons from — > ir Q ir° decays. The relation between 
A'£ — > 7T 7r 7r° and A'£ — > ir°ir° is similar to the rela- 
tion between A" — > 7r°7r° and A'£ — > iPvv. By taking 
their branching ratios into account, it was found that the 
low-mass region could to be reproduced after reducing 
the A° — > 7r 7r°7r° yield by several orders of magnitude. 
Reproduction could not be achieved without good sim- 
ulations of small signals from the detector, to which a 
tight veto was applied. This achievement indicates that 
a similar direction will be promising in the next step. 

The E391a experiment was also able to study other 
decay modes including 7r°'s in the final state, such as 
K\ -> n°n°vD and K° L -» t: tt X (X -» 77). The 
7r reconstruction method and the hermetic veto system 
are also valid in the analysis of these modes. Results 
based on portions of the data have already been pub- 
lished [H, Hi| , and further studies are in progress with 
the entire data set. 

The next step, the KOTO experiment [35[ at the new 
J-PARC accelerator (3(|, is now in preparation. Most 
of the improvements pointed out above are being imple- 
mented. The new beamline has been constructed; runs to 
evaluate properties of the beamline started in October, 
2009. 
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